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Abstract. We study energy relaxation in a phenomenological model for polymer built from rheological 
considerations: a one dimensional nonlinear lattice with dissipative couplings. These couplings are well 
known in polymer's community to be possibly responsible of /3-relaxation (as in a Burger's model). After 
thermalisation of this system, the extremities of the chain are put in contact with a zero-temperature reser- 
voir, showing the existence of surprising quasi-stationary states with non zero energy when the dissipative 
coupling is high. This strange behavior, due to long-lived nonlinear localized modes, induces stretched 
exponential laws. Furthermore, we observe a strong dependence on the waiting time t^ after the quench 
of the two-time intermediate correlation function C(tu,+t, t TO ). This function can be scaled onto a master 
curve, similar to the case of spin or Lennard- Jones glasses. 
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1 Introduction 

An important challenge in polymers physics nowadays is 
the understanding of non-equilibrium glassy state and as- 
sociated slow relaxations [[[]. At temperature below the 
liquid-to-glass transition temperature, the structural re- 
laxation time depends indeed on the time spent in the 
glassy phase, the so-called waiting time t w : this is the 
aging effect Low- frequency mechanical spectro- 

scopies show two different slow relaxations, namely a and 
(3, corresponding to two different temperatures T a > T^ 
which depend on the mechanical excitation frequency; T Q 
is shown to be very close to the liquid-to-glass transition 
temperature [Q. Both relaxation modes are characterized 
by a maximum of dissipation but are surprisingly differ- 
ently affected by physical aging: strongly for the first one 
but very weakly for the second. In spite of many studies, 
the microscopic origin of liquid-to-glass transition and of 
these relaxation modes are not yet well established f|. 
The a relaxation, associated to the liquid-to-glass tran- 
sition, involves cooperative motions inducing an increase 
of the viscosity, whereas localized motions would be re- 
sponsible of the [3 relaxation with a clamping of degrees 
of freedom (for instance macromolecule rotation) near the 
temperature Tp. 



Send offprint requests to: Franck Gobet 



Some phenomenological models (for instance Burger, 
Maxwell, Kelvin- Voight or Zener's models) have been de- 
veloped using linear viscoelastic theory and describe 
easily some aspects of relaxation in spite of their simplic- 
ity. These systems (see for example Burger's model in Fig. 
|l|a) correspond to a succession of springs and viscous ele- 
ments (dashpots) set in parallel or in series, to modelize 
interactions between and inside polymers. As we will see 
in Section 2, in the case of Burger's system, the dashpot 
with viscous parameter 71 (which modelized viscous inter- 
action between molecules) is linked to a relaxation (and 
therefore to liquid-to-glass transition), whereas the second 
piston (which modelized a dissipative interaction inside 
the molecules) induces (3 relaxation. Let us notice that 
these simple dynamical models don't exhibit any physical 
aging. One way to describe and explain this phenomenon 
could be that of taking into account nonlinear interactions 
inside macromolecules. 

Indeed, one of the most important feature of nonlinear 
systems is the localization of vibrational energy (named 
discrete breather) that modifies strongly the energy relax- 
ation P,^0[ and induces non-equilibrium dynamics. Brea- 
thers are known to have very interesting dynamical prop- 
erties and have been used in explaining a variety of physi- 
cal and biophysical phenomena. For instance in nonlinear 
systems where breathers become mobiles, they could con- 
tribute directly to the energy transfer and modify its re- 
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Fig. 1. Sketch of the Burger's linear model (panel a), and of 
our non-linear dissipative model for polymers (panel b). 



laxation properties in a nonexponcntial dependance. This 
interesting phenomenon has been invoked in several phys- 
ical settings as DNA molecules [p), hydrocarbon struc- 
tures flS^I , targeted energ y t ransfer between donors and 
acceptors in biomolecules [n3l . When the coupling is much 
smaller than the non linearity, the presence of essentially 
pinned long-lived breathers in nonlinear systems blocks 
the energy propagation The macroscopic manifesta- 
tion of this phenomenon is a very slow relaxation of the 
total energy, reminiscent of the long lifetime of metastable 
states in glassy systems observed after a quench. 

The case of relaxation phenomena in nonlinear lattices 
with bulk dissipation (as in glassy polymers) has received 
much less attention, though the experimental systems be- 
long to this class |l4| . Indeed, the (3 relaxation in glassy 
polymers, which is characterized by a maximum of dis- 
sipation, may be induced by the clamping of degrees of 
freedom. In spite of this dissipative effect, physical aging 
of this system is very slow. Understand this strange be- 
havior is one of the aims of this paper. To accomplish this 
goal, we will examine the influence in energy relaxation of 
a dissipative coupling by performing numerical studies on 
a phenomenological model of polymers. The system, de- 
scribed in Section 3 and pictured in Fig. ^b, corresponds 
to particles coupled via elastic and dissipative interactions 
(as in Burger's model for the molecular modelization); in 
addition, each particle is submitted to an on-site nonlin- 
ear potential to take into account interactions between 
different polymers including several contributions such as 
hydrogen bonds linking and repulsion of two chains (short- 
range steric restrictions). For a fixed viscous parameter 7, 
we examine the relaxation of energy when, after thermal- 



ization, the ends of the chain are placed in contact with 
a zero-temperature reservoir. Results show different kinds 
of energy relaxation regime which depend strongly of the 
dissipative terms: in particular, we show that the system 
can relax very slowly in spite of high dissipative couplings! 

2 Viscoelastic systems and a, (3 relaxations 

In many experimental studies, isochronal dynamical me- 
chanical spectrometry has been performed at different fre- 
quencies Q] . With these measures, it is possible to map the 
loci of a and relaxations (defined by maximum of the 
imaginary parts of shear modulus G=G'+iG" ) versus tem- 
perature and frequency. Typically, the a-relaxation pro- 
cess exhibits a clearly non-Arrhenius behavior above the 
liquid-to-glass transition T s , whereas an Arrhenius behav- 
ior is observed for temperature below T g , when the system 
is out of thermodynamical equilibrium. The second pro- 
cess (3 appears at higher frequencies than a-process with 
also an Arrhenius behavior. 

It is well-known in polymer's community that some as- 
pects of the rheology of glassy polymers can be described 
easily by usingphenomenological models from viscoelastic 
linear theory [pl|7l||. It can be assumed that the defor- 
mation of the polymer is divided into elastic and viscous 
components and can be described by a combination of 
Hooke and Newton's laws. These systems correspond to 
springs and dashpots, either in series or in parallel. The 
Burger's model (one element is described in Fig. |l|a) seems 
to be a good candidate to illustrate the rheological behav- 
ior of polymers in glassy state. One element is made of two 
springs in series to modelizc the chain of a macromolecule. 
In addition, one part of this chain is coupled in paral- 
lel with a dashpot to take into account the clamping of 
some degrees of freedom at low temperature, inducing the 
/3-relaxation. The second dashpot in series with elastics 
components illustrates the viscous behavior of the system 
when it becomes glassy and induces the a-relaxation. 

In this model, the shear modulus G=G'+iG" is mea- 
sured by applying a sinusoidal force, with frequency u>, 
at one end of a chain with 100 elements and by fixing 
the other extremity. G is the ratio between the applied 
force and the relative displacement of the chain extremity 
obtained in the stationary limit. The shear modulus is de- 
termined in a temperature window by using two Arrhenius 
laws for the viscosity 71 and 72: 



and 72 = 10~ 5 7i. Furthermore, as in all viscoelastic linear 
studies (see for instance references ^,0,||)) there is no in- 
erty inside the system: each connection between springs 
and dashpots is considered as a particule with a mass 
equal to zero. 

In Fig. we report the evolution with temperature of 
the shear moduli G', G" and loss tangent tan0=G"/G' 
for parameters given in the corresponding caption. These 
curves show two distinct regions of rheology, namely for 
increasing temperature, the /3-process (weak decreasing 
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of G' and, maxima of G" and tan0) and the cv-process 
(strong decreasing of G' and maximum of G"). This be- 
havior is qualitatively similar to those observed experi- 
mentally ||. However, some differences appear. The loss 
tangent exhibits an increase with temperature for the a- 
process, but we don't observe a maximum unlike in ex- 
perimental measurements. Furthermore, for temperature 
higher than the liquid-to-glass transition temperature, a 
non-Arrhenius behavior is observed experimentally for the 
cv-relaxation. This characteristic is not take into account 
in the evolution of 71 . 
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Fig. 2. Values of the real (G') and imaginary (G") parts of 
shear modulus and of tan0; they are plotted versus tempera- 
ture for the Burger's model of Fig. |l|a. The shear moduli are 
divided by the modulus limit of low temperature Go- Parame- 
ters used in the calculations arc k=0.1, 7?=10 6 , E/k 6 =9.7 and 

^=icr 2 . 



In spite of these differences, the major overriding con- 
clusion for this section is the possibility to observe rheo- 
logical characteristics, like /3-process, of the glassy state by 
considering a viscous coupling inside the molecule (linked 
to the rotation of terminal groups or other side chain mo- 
tions ||). This modelization is however very concise and it 
is not possible to describe any aging phenomena because 
of the linearity of components. 

In the next part, we will describe a phenomenologi- 
cal model including dissipative couplings and nonlinear 
interactions that may modelize more precisely one chain 
of a macromolecule. The main question we would like 
to address is the possibility to observe a long lived non- 
equilibrium state in nonlinear systems despite dissipative 
terms. In order to solve this problem, we simplify the sys- 
tem: the initial thermalized nonlinear system is put in a 
zero-temperature bath with a fixed viscous parameter in- 
dependent of temperature. 



3 The phenomenological dissipative nonlinear 
model 

We consider a one dimensional chain of N=200 anhar- 
monic oscillators, with a nonlinear on-site potential V(x), 
with free ends and nearest-neighbor elastic coupling po- 
tential (the coupling being k). The sketch of the chain 
is reported in Fig. lb. In analogy with Burger's model 
and in order to modelize the clamping degree of freedom, 
we consider for one half of nearest-neighbors a dissipa- 
tive coupling (7 is the dissipative parameter) in parallel 
with clastic coupling. For the on-site potential V(x), de- 
scribing the interactions between two molecules (hydrogen 
long-range attraction and steric short-range repulsion) , we 
have chosen the Morse potential 



V(x) = i[l 



-XI 2 



(2) 



which has the appropriate shape to describe the strong 
repulsion when the chains are pushed toward each other 
(x<0) and the vanishing interaction when the chains are 
pulled very far apart (x^>l). Each end oscillators of our 
chain can be also submitted to an additional damping 
force. The equations of motion of this chain are given in 
dimensionless form by: 



X 2 i 



and 



dV 
dx 2 i 
-7(x 2 i 



dV 



X 2 i+1 = 



k(2x 2i - x 2i+ i - x 2i -i) 
x 2i _ 1 )-7'x 2i ,5(2i- 200) 

k(2x 2i+ i - x 2i - x 2i+2 ) 



(3) 



<9x 2i+ i 

~7(x 2 i+i - x 2i+2 ) - 7'x 2i+1 (5(i) (4) 



where x„ is the dimensionless displacement of the nth os- 
cillator from equilibrium, x n its velocity and 8 the Kro- 
nccker delta function. The mass of the oscillators is set 
to unity by appropriately renormalizing time units. The 
equations of motion have been integrated using a fourth 
order Runge-Kutta method. 

To study energy relaxation, we consider k=0.01 and 
we initially thermalize the system at temperature T=l by 
using Nose-Hoover thermostats |fT6f . This temperature is 
much higher than the critical temperature (T c =0.2) of the 
"order-disorder" transition which characterizes the non 
dissipative model 7=0 (for more detail see references [ 17| 
|l8t). Then, in average, the kinetic energy per site is higher 
than the depth of Morse potential (which is equal to 0.5 
in our arbitrary units) and we can consider our system 
as in an initial "liquid" state: all the particules are in the 
plateau of the Morse potential. In other words, there are 
no interaction between chains of molecules. The thermal- 
ization procedure is performed with 7=7'=0 by using a 
chain of three thermostats to provide a good exploration 
of the phase space |l6| . 

After thermalization, the connection with the heat bath 
is turned off and the lattice is connected to a zero temper- 
ature reservoir via the damping term with 7'=0.1. How- 
ever, we explore the interval [0,100] for the dissipative 
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parameter 7. As our purpose is to search for long lived 
non-equilibrium state in spite of dissipative couplings, we 
prefer to simplify the problem and not to consider a tem- 
perature dependence of 7 like in the Burger's model of 
section 2. The moment of connection with the zero tem- 
perature reservoir is chosen as the origin of time. At each 
step of the integration of equations (2) and (3), we evalu- 
ate the total lattice energy: 



21 Ml 



ix 2 
2 1 



V(x- 



1!)!) 



(5) 



and consider the symmetrized local energy per site: 

Ei = ix2+V(x i ) + ^(x i _ 1 -x i ) 2 + ^(x i -x i+1 ) 2 . 

(6) 

Total energy is expected to decrease with time and con- 
verge to a zero value of the "frozen" state at equilibrium. 



4 Relaxation of the thermalized system 

In Fig. ||, we report the total lattice energy divided by the 
initial energy versus time for various viscous parameter 
7. For 7=0, a clear non exponential decreasing energy is 
observed, consistent with Tsironis and Aubry's results Q. 
This long-tail relaxation behavior was shown by these au- 
thors to be connected to the presence of long-lived non 
linear localized modes that are relatively mobile. If we 
consider a small dissipative coupling (7=10 -3 ) inside the 
lattice, we see that total energy decreases faster than pre- 
viously. This small dissipative coupling induces a strong 
modification of energy relaxation to the "frozen" state 
(normalized energy closed to 0) for a time smaller than 
10 4 (in arbitrary units). This emphasizes that the dissipa- 
tive couplings change strongly the relaxation mechanisms 
and can induce a fast relaxation regime. 

For 7 higher than 0.1, a new surprising feature is ob- 
served: normalized energy seems to be blocked with a very 
slow decrease for long time whereas the dissipative pa- 
rameter is higher! The system seems to evolve in a quasi- 
stationary state that is neither a "frozen" state (the nor- 
malized energy is clearly different from 0) nor a "liquid" 
state" (energy is too low). We can also notice that, at a 
given time, the energy of this quasi-stationary state in- 
creases with 7. This situation seems to be qualitatively 
similar to polymer systems where a very long-lived non- 
equilibrium glassy state is observed with very slow phys- 
ical aging for temperature smaller than T Q (the viscous 
parameter which follows an Arrhenius law is higher for 
lower temperature). 

In previous relaxation studies in non linear lattices 
where blocking energy was observed, it has been shown 
that such behavior may be induced by long-lived breathers 
||. In order to examine more precisely the present situa- 
tion, we report in Fig. ^ the spatiotemporal energy land- 
scape of the lattice by plotting the local energy E, in each 
lattice site for 7 = 10 _3 ,10 _1 and 10. Time advances 
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Fig. 3. Normalized total energy decay as function of time for 
various viscous parameter 7. Initially each array is in thermal 
equilibrium at the temperature T=l. 



along the y axis until t=10 4 and a gray scale is used to rep- 
resent the local energy with darker shading corresponding 
to more energetic regions. For 7=10 -3 , we can see two 
kinds of energy relaxation: on the one hand, there is a dis- 
sipation of the energy inside the lattice, characterized by a 
" fibrous structure" of the local energy landscape in Fig. [|c 
and, on the other hand, we observe a dissipation of mo- 
bile breather via surface damping characterized by "dark 
oblique lines" . For 7 >0.1, we notice the clear presence of 
pinned long-lived breathers, responsible of the energy re- 
laxation blocking. This localisation of energy is observed 
after a short time where energy not only decreases via 
the surface damping but also via a dissipation inside the 
lattice (see for example the landscape for 7 = 10 _1 and 
t<2500). Furthermore, these states have a very long lived 
time and are still observed for t higher than 10 5 . 

We have mentioned previously that, at early times, 
phonons and mobile breathers dissipation takes place be- 
fore pinned breathers relaxation. Typically the hierarchy 
of relaxations processes may be classify in a sequence 
of characteristic times tj <t?> <... (ljj, where the en- 
ergy relaxation corresponds approximately to exponential 

Let us introduce 



or a stretched exponential decay |20 



the phonon relaxation time t* defined by E(t*)/E(0)=0.5 
which is characteristic of the phonon and mobile breather 
dissipation inside the lattice. We have reported in Fig. || 
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Fig. 4. Evolution of the local energy Ej along the chain for 
various viscous parameter 7=10 (panel a), 0.1 (panel b) and 
0.001 (panel c). The grey scale goes from Ei=0 (white) to the 
maximum Ei-value (black). 
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Fig. 5. Phonon relaxation time t*, defined by E(t*)/E(0)=0.5, 
for various values of the viscous parameter 7. 



is close to 0.5 at time t=2.10 4 , and the associated sites are 
not linked with the on-site Morse potential. This quasi- 
stationary state is clearly not a frozen state, since some 
parts of a molecular chain are "hot" but do not interact 
with other ones. 



the evolution of this relaxation time versus the viscous pa- 
rameter 7. We clearly see, in this first step, that energy 
decreases faster for 7 between 0.1 and 1. In the case of 
a linear lattice, the maximum of dissipation is predicted 
to correspond to a value of 7 which verifies jlu ~ K + k 
where K corresponds to the coupling constant of the lin- 
earized Morse potential (K=l in our case) and uj is the fre- 
quency of the phonon band that starts at the frequency 
VK=1 and extends to \/K + 4k~1.02. Therefore, in the 
case of linearized oscillator, we expect a maximum of dis- 
sipation for 7 close to 1 in agreement with what is re- 
ported in Fig. |^. It is thus very surprising to observe quasi- 
stationary states in a second regime whereas dissipative 
effects are important. We have to notice that this para- 
doxal effect is very similar to the situation observed in the 
case of P relaxation of polymers in glassy state: a very slow 
physical aging in spite of the maximum of dissipation Q . 

In Fig. ||, we report the local energy in the lattice for 
various parameter 7 and at a given time t=2.10 4 much 
higher than the characteristic time of phonons and mobile 
breathers dissipation. For 7<0.1, the local energy per site 
is close to as seen previously in Fig. |[ the system with 
small dissipative coupling reached thus rapidly an equi- 
librium frozen state. For 7>0.1, we clearly see the local- 
ization of energy, as long-lived breathers, in two nearest- 
neighbors oscillators dispatched in some areas of the lat- 
tice. We have verified that the two considered nearest- 
neighbors are coupled via piston: when 7 is high enough 
this coupling is clamped. In fact, the phase displacement 
of velocities of the both particles is very small, inducing 
a very slow dissipation and then referred a very long non- 
equilibrium state. The energy of these long-lived breathers 
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Fig. 6. Instantaneous local energy Ei along the chain at time 
t=2.10 4 for various dissipative parameter 7. The dashed lines 
show the limit of Morse potential for large displacement x^. 



For t^>t* and 7>0.1, the energy relaxation of the quasi- 
stationary state can be very well fitted by a Kohlrausch- 
Williams- Watts function or stretched exponential func- 
tion: 



E(t) = E b (0)e 



-(at) b 



(7) 
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where the coefficients a, b and Eb(0) are 7 dependent. 
Eb(0) can be qualified as the total energy of pinned brea- 
thers at time t=0. The Kohlrausch |2l[] exponent b is a 
parameter measuring the deviation from a single exponen- 
tial form (0<b<l). In Fig. we have reported the evolu- 
tion of [-lnE(t)/Eb(0)] versus t with logarithmic scales for 
various parameter 7>0.1. We can see a very clear linear 
dependance for high t (after phonon and mobile breather 
dissipation), which attests that this long-lived breather 
energy relaxes as a stretched exponential. 
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Fig. 7. Log[-ln(E(t)/Ef,(t))] versus log t for various dissipative 
parameter 7. The dashed lines are the stretched exponential 
fits. 



We see also in these figures that the slope of the straight 
lines depends on the viscous parameter 7. In order to ex- 
amine more precisely this dependence, we report in Fig. || 
the Kohlrausch exponent b versus the viscous parame- 
ter 7. We clearly see a maximum of the exponent b for 
values of 7 close to 0.5. The b-value is equal to 0.82 and 
then the energy relaxation differs from a pure exponen- 
tial decay. Furthermore, this figure shows that the pinned 
breathers relaxation is slower for higher parameter 7. It is 
therefore qualitatively similar to the glassy state where 
a slower aging phenomena is observed for slower tem- 
perature (and then higher viscosity as we can see in the 
Burger's model of section 2) 0. 
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Fig. 8. Kohlrausch exponent b for various viscous parameter 7. 



5 Out of equilibrium dynamic correlations 

Studies of noncquilibriu m sy stems like spin, structural or 
Lennard- Jones glasses |5, 22 have shown that the nonequi- 
librium dynamics of the previously described states could 
be much efficiently characterized by two-time correlation 
functions of the form: 



C(t w + t,t w ) = (A(t w - 
-(A(t v 



t)A(tw)> 
+ t)}(A(t w )) 



(8) 



where A is a microscopic observable, and t w is the "wait- 
ing time" i.e., the time elapsed after the quench. Brack- 
ets denote an average over different initial configurations 
at temperature T. At equilibrium, this two-time quantity 
satisfies time translation invariance and then depends only 
on the time t. On the other hand, in out of equilibrium 
situations, such equilibrium property is not verified: this 
function depends on the waiting time t w ("aging effect"). 
The correlations functions for large times are expected to 
scale in the form: 



C(t w + t, t w ) = CsT(t) + Cag 



£(tw + t) 



(9) 



The first term describes short time dynamics that does not 
depend on t w and has the equilibrium form. The second 
term, or aging part, depends only on the ratio £(t w +t)/ 
£(tu>) where £(t) is a monotonous increasing function of t. 
In a lot of cases £(t)oct or oct" so that the aging part is 
simply a function of t/t w and exhibits a master curve (see 
for instance, experiments on thermorema ncnt m agnetiza- 
gels 



5,26|) 



tion ]23|, gels ]2J| or particle suspensions 

In this study we have considered the microscopic ob- 
servable A(t)=(2i=i Xi(t))/200 that is the mean deforma- 
tion per site of the chain at time t. Numerical calculations 
have been done in the case of a quench of the system 
with a viscous parameter 7=10 from temperature T=l 
to T=0. Two-time correlation functions are obtained for 
various waiting time by considering 1 1 different initial con- 
figurations. Furthermore, in order to make a quantitative 
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comparison, we prefer to calculate correlation normalized 
by C(t„,0. 

In Fig. |9|, we have reported the evolution of normalized 
two-time correlations function C(t l0 H-t,t UJ ) versus time t 
for different waiting time (we consider t^^t* in order to 
study only long lived non equilibrium states after phonons 
and mobile breathers dissipation). The behavior of C clearly 
emphasizes the lost of time-translation invariance and the 
dependence on the waiting time t w . This figure also shows 
that the dynamics can be decomposed into two time scales: 

(i) at short time separation (t<20) correlation function 
doesn't depend of t w and is equal to the value expected 
at equilibrium (C=l at T=0). 

(ii) the decay from this value toward zero arises in 
a second time scale that clearly depends on t w : the sys- 
tem does'nt reach equilibrium within the time window ex- 
plored in the simulation. Furthermore we can notice that 
the larger the waiting time, the longer it takes the system 
to forget the configuration at time t^. This behavior is 
typical of aging effect |1 . 
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Fig. 9. Normalized non-equilibrium correlation function 
C(t t0 +t,tti,) for fixed value t w versus time t. The viscous pa- 
rameter is 7=10 and from bottom to top the waiting times are 
t u -5844, 9300, 14700, 23137, 36321, 56920, 71225. 

Guided by equation (^) we test the scaling assump- 
tion for long times in Fig. [l^, where normalized correla- 
tion functions are reported versus normalized time t/t w . 
The different curves can be superimposed, indicating the 
validity of the scaling ansatz and the existence of a mas- 
ter curve. This striking feature is observed in many non- 
equilibrium systems p3|[2^ ,p5|,p6| and is very similar with 
comparable studies on spin glasses || or Lennard-Jones 



glasses |22| . The physical origin of this universal t /t w scal- 
ing is, at this day, an open question. Kob and Barrat sug- 
gest in reference |22j that it could be induced by a simi- 
larity of the geometry of phase space of these systems in 
spite of differences in microscopic dynamics. 

Finally, we would like to point out that Fig. || shows 
also the violation of the dissipation-fluctuation theorem 
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Fig. 10. The data of Fig (J versus t/t w : master curve 



(FDT) . It seems to be a characteristic of non-equilibrium 
system as observed numerically in domain growth process 
[ p7[ , Lennard-Jones glasses pq| or slow granular rheol- 
ogy p9| and experimentally for dielectric measurements 
in colloidal glasses [^0[ , supercooled fluid |?l| and polymer 
glasses |52||3j|. Let us consider the response R(t u) -|-t,t t0 ) 
to a field h conjugated to observable A. At equilibrium, 
response to this external field is given by FDT i.e., at 
equilibrium, the following equation is verified: 



R(t w + 1, t v 



<5(A(t w + t)> 



h=0 



and then: 



6>C(t w + t,t v 

at 



<5h(t w ) 
1 <9C(t w + t,t w ) 
T dt 



— — T R(t w + 1, t w ) 



(10) 



(11) 



where T is the bath temperature. 

If we consider a quench at a zero-temperature bath (as 
in our case), then for satisfying FDT, normalized correla- 
tion function has to be constant: 

C(t w + t,t w ) = C(t w ,t w ) = l , Vt>0. (12) 

Evolutions of correlation functions reported in Fig. |^ show 
clearly that C(t UI +t,t u; ) is not a constant with time t: 
fluctuation-dissipation theorem is violated in this non- 
equilibrium system as seen in other glasses ]27|,p8|, p9| , |3C| , 
[jf| . [32||33j . Moreover, the larger the waiting time, the longer 
it takes the system to violate the FDT. 



6 Conclusion 

Using simple rheological considerations and results from 
Burger's model, we built a phenomenological model of a 
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polymer-chain. This system is a one dimensional nonlinear 
lattice characterized by dissipative couplings. The energy 
relaxation studies of this thcrmalizcd system show that for 
sufficiently large viscous parameter 7, it is possible to ob- 
serve nonequilibrium quasi-stationary states in spite of the 
short characteristic time of phonon dissipation! This sur- 
prising behavior is due to a chain auto-organisation which 
minimizes energy dissipation, inducing the clamping of 
some degrees of freedom and forming long-lived pinned 
breathers. 

Moreover, this very slow energy relaxation can be fit- 
ted by stretched exponential laws, ubiquitous in glassy 
polymer aging properties. Another similarity with these 
physical systems is that this aging phenomenon is slower 
when the viscous parameter 7 is higher. Furthermore, the 
two-time correlation function C(t w +t,t w ) shows a strong 
dependence on the waiting time and can be scaled onto a 
master curve by considering the evolution versus normal- 
ized time t/t„ ; . 

By these aspects (clamping of degree of freedom, long- 
lived nonequilibrium state, stretched exponential decays, 
violation of time-translation invariance and master curve), 
the glassy nature of this simple dissipative non linear lat- 
tice is very similar to those observed in glassy polymers. 
Beside its interest for nonlinear physics, this model is 
presumably an alternative to study complex systems like 
glassy state polymers: we have now to push our inves- 
tigations further by examining other properties like de- 
pendences on cycling temperature, evolution with wait- 
ing time of the shear moduli. Work along this line is in 
progress. 
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